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We study the spectrum of fermions hopping on a chain with a weak incommensuration 
close to dimerization; both q, the deviation of the wave number from ir, and 5, the 
strength of the incommensuration, are small. For free fermions, we use a continuum 
Dirac theory to show that there are an infinite number of bands which meet at zero 
energy as q approaches zero. In the limit that the ratio q/S — > 0, the number of states 
lying inside the q = gap is nonzero and equal to 25/% 2 . Thus the limit q — » differs from 
q — 0; this can be seen clearly in the behavior of the specific heat at low temperature. 
For interacting fermions or the XX Z spin-1/2 chain, we use bosonization to argue that 
similar results hold. 



PACS number: 71.10.Fd, 75.10.Jm 



One-dimensional lattice models with incommensurate hopping elements or on-site potentials have been 
studied for many years from different points of viewETH. Many unusual properties of the quantum spectra 
and wave functions have been discovered for various kinds of aperipdicityu. Physically, such models have 
applications to several problems such as incommensurate crystalsO and semiconductor heterojunctionsEl. 
If the incommensuration is close to dimerization (wave number tt), then the models have an additional 
interest in the context of metal-insulator transitions or spin-Peierls compoundsLrtl. In this paper, we 
study a model of fermions with a hopping which has a small incommensurate term close to dimerization. 
We will show that this has some unusual properties, particularly inside the gap which exists exactly at 
dimerization. The major surprise is that states appear inside the gap as soon as we move away from 
dimerization. These states carry a finite weight, and therefore contribute to quantities like the specific 
heat at low temperatures. The bulk of our discussion is for free fermions or the XY spin-1/2 chain for 
which quantitatively accurate results can be obtained by analytical methods. However, we will argue at 
the end that the same results should hold for interacting fermions or the XX Z spin-1/2 chain, which is 
the more general and interesting case. 

We begin with the following Hamiltonian for noninteracting and spinless fermions on a lattice 



H = 

2 



2 ^ n ( C ™ Cn + 1 C n+l C n ) A* c li c n i 

n n 

J n = 1 + 5 cos (n + q)n , (1) 
where we will assume that <5 <gC 1 and q <C n. We set the chemical potential \x = as we are interested in 
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energies close to zero. If 5 = 0, the model can be easily solved; the dispersion relation is E(k) = — cosfc. 
In the ground state, all the states with momenta lying in the range [— kp, kp] are filled, where the Fermi 
momentum kp = ir/2. The Fermi velocity is equal to sin kp = 1. 

If 5 but 9 = 0, i.e., with dimerization, the model can still be solved analytically. There is an energy 
gap extending from —5 to S. Let us now consider nonzero values of q. [Since cos(jr-\-q)n — cos(qn) cos(7m), 
we see that a small value of q is equivalent to dimerization term whose amplitude 5 cos(qn) varies slowly 
with n]. For any rational value of q/ir = M/N, where M and N are relatively prime integers, we have a 
periodic system with period P equal to N if M + N is even and 27V is M + N is odd. The one-particle 
spectrum of Eq. ([IJjpCan be found by computing a 2 x 2 transfer matrix M(q, E) obtained by multiplying 
together P matricesa, 



If |tr M(q, E)\ < 2, the energy E is allowed in the spectrum; otherwise that energy is not allowed. By 
following this method and sweeping through a large number of values of q and E, we obtain the picture 
of energy bands and gaps (shaded and unshaded regions, respectively) shown in Figs. 1 and 2, taking 
5 = 0.05. (We have sampled the energy values with a resolution of dE = 10~ 6 ). We have scaled both E 
and q in units of 5 because the pictures turn out to depend only on those two ratios. We immediately see 
that the pictures look much more complicated than the situation in which q is exactly equal to zero; in 
that case, there is precisely one gap extending from E/S = to 1. [We will consider only positive values 
of q and E since it can be shown that the spectrum is invariant under either q — > — q or E — > —E. Thus 
the energy bands shown in Figs. 1 and 2 should be understood as continuing to negative values of both 
E and q by reflection about those two axis] . 

Fig. 1 shows that as q/S increases, the gaps shrink rapidly. We will show that this can be understood 
using perturbation theory to n th order, where n is an odd integer, and that the gaps shrink as (6/q) n . We 
have therefore labeled the gaps in Fig. 1 by the integers n=l,3,5,.... More interestingly, we observe that 
all the bands approach the origin (q, E) = (0, 0). The widths of the low-lying bands vanish exponentially 
fast as q/5 — ► as we will argue from a WKB analysis. This makes it impossible to find these portions of 
the bands by looking for energies satisfying |tr M(q,E)\ < 2 using any reasonable energy resolution dE. 
We therefore find these thin portions of the bands by looking for minima of trM 2 (g, E) as functions of E; 
these minima yield single points in energy which lie within a distance of dE of the band. We are thus able 
to find the locations of the thin regions without having to use a resolution smaller than dE = 10~ 6 . In 
Fig. 2, we show these points (thin regions) which smoothly join on to the wider regions of the bands. We 
also show six curves which are the analytical results of a low-energy theory of the model which involves 
solving a 1 + 1-dimensional Dirac equation in a periodic potential. We will also prove that the number 
of bands in the region < E/S < 1 increases as 26/irq as q/S — ► 0, while the number of states in each 
band is equal to q/n when normalized appropriately in the thermodynamic limit, i.e., the number of sites 
L — > oo. This will show that the number of states lying inside the q — gap is finite and equal to 25 /ir 2 
in the limit q —> 0; this implies that q = is a rather singular point. 

To begin with, let us understand the gaps for large values of q/5 by using perturbation theory about 
5 = 0. If q and 5 are both small, the states close to zero energy are dominated by momenta near ±7r/2. 
Since the incommensuration has Fourier components at momenta ±(ir + q), we see that the gaps above 
zero energy result from the breaking of the energy degeneracy between the two states with momenta 
k\ = 7r/2 + nq/2 and k n+ \ = —tt/2 — nq/2, where n = 1,3,5, .... These two states are connected to 
each other through the n — 1 successive intermediate states with momenta = ~ir/2 + (n — 2)q/2, 
kj, = ir/2 + (n — 4)g/2, k^ = —ir/2 + (n — Q)q/2, k n = tt/2 — (n — 2)q/2, since we must have fc^+i — k{ = 
±(7r + q) mod 2tt. Since all the relevant matrix elements are equal to (5/2, and the energy denominators 
are of the order of q each, we see that the two-fold degeneracy (at the energy E\ = E n+ i ~ nq/2) gets 
broken at the n th order in perturbation theory to produce a gap of the order of 5 n /q n ~ 1 . To be explicit, 
we find that the upper and lower edges of the two gaps labeled by n = 1 and 3 in Fig. 1 are given by 




p 



(2) 



n=l 
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E±(l) q_ 1 

S 25 2 ' 

g±(3) 3g jW 

<5 25 + 16g 32 g 2 ' [ ' 

and we have verified that this agrees well with the numerically obtained boundaries. We will show 
elsewhere that the general formula for the gap AE(n) — E + (n) — E-(n) labeled by the odd integer 
n = 2p + 1 is given by 

Ag(n) (S/Ag)^ 

5 (pi) 2 { 1 

to lowest order in 5. 

We can count the number of states in the band lying between gaps n and n + 2 as follows. For 5 = 0, let 
us normalize the number of states so that the total number is 1; since the momentum k goes from —tt to tt 
in that case, the density of states in momentum space is 1/27T. If we now turn on a very small value of 5, we 
find that the band lying between the gaps n and n + 2 is made up of linear combinations of the states lying 
between the two momenta intervals [7r/2 + nq/2, ir/2 + (n + 2)q/2] and [— tt/2 — nq/2, —ir/2 — (n + 2)q/2]. 
The total number of states in these two intervals is 2(g/2-7r) = q/-n. Thus each band contains q/it states. 
Now, this number cannot change if we change 8, and we can therefore use the same number below in the 
opposite limit where q/5 is small. 

Let us examine the more interesting regime where q/8 is small and E/5 < 1. For analyzing this, it is 
useful to consider the continuum limit. In this limit, the Fermi field ^(x) = c n can be written as 

7T 7T 

t/j(x) = ipR(x) exp (i-^x) + iPl{x) exp (-i-x) , (5) 



where ipR and ipL denote the right- and left-moving fields, respectively; they vary slowly on the scale 
of the lattice spacing (which is set equal to 1 here). We substitute (||) in Eq. ([]]) and drop terms like 
exp(±i7ra;) which vary rapidly. We then find the following Dirac-like equations for the two time-dependent 
(Heisenberg) fields 

i ( d/dt + d/dx ) — iScos(qx) tpL = , 

i ( d/dt — d/dx ) ipL + iScos(qx) tpR = . (6) 

Eq. (|^) can be solved by defining ip± = ipR ± ipL which satisfy the second-order equations 
[d 2 /dt 2 - d 2 /dx 2 + 5 2 cos 2 (qx) - 6qmn(qx) ] ip+ = 0, 

[d 2 /dt 2 - d 2 /dx 2 + S 2 cos 2 (qx) + 6qsm(qx) ] V- =0. (7) 

It is sufficient to solve one of these equations, say, for ip + , since ip- is related to The energy spectrum 
can be found by solving the time-independent equation 

[ - d 2 /dx 2 + 6 2 cos 2 (qx) - 5qsm(qx) ] tp + = E 2 ip + . (8) 



Eq. (||) has the form of a Schrodinger equation (with "energy" E 2 ) m the presence of a periodic 
potential. The potential is similar but not identical to Mathieu's equationB. We therefore again expect 
bands to form. To simplify the notation, let us shift x by ir/2q and then scale it by a factor of q to make 
the period equal to 2it. We then get 

S 2 S E 2 

[ — d 2 /dx 2 + —sin 2 a; — -cosx] ip+ = — -0+ . (9) 
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By Floquet's theorems, the solutions must satisfy tp+(x + 2ir) = e l9 ip + (x), where 9 goes from to 
±7r from the bottom of a band to the top. We observe that there is an exact zero energy state with 
i/>+(x) — exp( | cos x ) and ?p-(x) — for which = 0. In general, the nonzero energy states can only 
be found numerically. But if q/S is large, the positions of the low- lying bands (with E 2 /q 2 <C S/q) can 
be found analytically as follows. For energies much lower than the maxima of the potential, we can 
ignore tunneling between the different wells to begin with. Near the bottom of a single well, we have a 
simple harmonic potential with small anharmonic corrections. Ignoring the anharmonic terms, we find 
the energy levels to be simply given by 

El = 2nqS , (10) 



where n = 0, 1,2, .... Next, we include the anharmonic corrections perturbatively. Up to second-order in 
perturbation theory, we get the more accurate expression 



E n = y/2^qd [ 1 



n q 
86 



5n 2 + 2q 2 
128 S 2 



(11) 



We can now include tunneling between wells; a WKB analysis shows that each of the above energies 
should split into a band whose width AE n is of the order of 6 times exp(— 28 / q). This explains why the 
bands rapidly become thin and shrink to isolated points as q — > in Fig. 2. In that figure, we have shown 
the six curves corresponding to n = 1 to 6 in Eq. (|i"l|); for n = 0, we simply get a straight line lying at 
zero energy. We see that all the curves agree extremely well with the numerical data, even up to E n of 
the order of 5 where the harmonic approximation breaks down and the band widths become noticeable. 

We will now prove that the number of states lying below the line E/5 = 1 is proportional to S. For 
large 5/q, we can use a semiclassical phase space argument to count how many states lie below any 
given energy. The "Hamiltonian" on the left hand side of Eq. (||) has the form p 2 + V(x), where 
V(x) = S 2 cos 2 (qx) — 5q sm(qx) . Hence the number of states up to energy E is given by the phase space 
integral 

Iff dxdp n/ 2 2 



v{E) = Z J J ^f @{E P V{X) ] ' (12) 

where we have divided by the length L of the system for normalization; the 0-function is defined to be 
1 and if its argument is positive and negative, respectively. On doing the momentum integral in (|l^), 
we get 

v{E) = ± dx - s/E 2 - V(x) 9( E 2 - V(x) ) , (13) 

2vr Jo 7r 

where we have used the fact that V(x) has period 27r/g. For small values of q/S, we then see that the 
number of states lying between E — and S is 

v{S) = % , (14) 



which is independent of q. Incidentally, this also yields the number of bands lying below E = S. Since 
each band contains q/ir states, Eq. ( pT[ ) gives the number of bands to be 25/wq. We have verified that 
this estimate agrees very well with our numerical results. 

The difference in the spectrum for q = and q — > should show up most clearly in the specific heat 
at temperatures T <C 6. (We set the Boltzmann constant equal to 1 for convenience). For q = 0, the 
energy gap from —d to +6 means that the specific heat vanishes exponentially at low temperature. But 
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for q — ► 0, the semiclassical expression in ([l3|) shows that the density of states at any energy E -C S is of 
the order of p(E) ~ E 1 /^. This implies that the low-temperature specific heat goes as T 2 /5 which is very 
different from the exponential behavior at q — 0. 

Let us pause to ask: why does the model in Eq. (|l|) show such different behaviors in the limit q — > 
and at q — 0? To answer this question, it is useful to generalize the form of the incommensurate hopping 
in (Q) from <5cos(-7r + q)n to <5cos[(7r + q)n + 77], where r\ is a phase lying between and 2ir. If q/ir is 
irrational (this is the generic case to consider if q is nonzero), the spectrum does not depend on r\. [If 
q/n = M/N is rational, the spectrum does depend on 77, but it varies less and less with r\ as N — > 00. 
In our numerical computations, we introduced a random phase 77 for each value of q and found that this 
had no noticeable effect on the two figures] . This insensitivity to 77 can also be seen from the continuum 
theory in Eq. (^J) where a potential of the form cos(qx + 77) can be transformed to cos(qx) by shifting x 
appropriately. However, exactly at q = 0, the hopping has a term like 8 cos rj cos(7rn) and the spectrum 
depends significantly on 77; for instance, there is a gap from E = up to | <5 cos 17 1 , and the number of 
states up to an energy E is given by 



v 



(E,ri) = - yjE 2 - S 2 cos 2 77 (15) 



7T 



for S cos 77 1 <£<gl. We may now consider taking an average of this number over all possible values of 
77. Thus, the number of states up energy E is given by Eq. (|l^) to be 

v{E) — C* ^E 2 - 8 2 cos 2 77 6( E 2 - S 2 cos 2 77 ) . (16) 
Jo 2tt 2 

This agrees with ( |T3| ) if we substitute 77 = qx. Thus, the limit q — > agrees with the point g = provided 
that we average over 77 in the latter case. However, this statement may not have any physical significance 
since an experimental system sitting at q — (the dimerized point) is more likely to choose a particular 
value of 77, rather than average over many values of 77. We will therefore choose 77 = henceforth. 

We will now argue that the difference between small q and q = (in particular, the presence of 
states within the q = gap) persists for the more interesting case of interacting fermions, i.e., for 
Luttinger liquids. To this end, consider adding a four-fcrmion interaction term like c^CnC^-yCn+i to 
the Hamiltonian ([!]). Equivalently, we can consider an XX Z spin- 1/2 chain governed by 



H = \ [ ( 1 + ^ C0S + ) ( S n S n+l + S l S l+l ) + A S n S n+l 



(17) 



which is related to the interacting fermion theory by a Jordan- Wigner transformation^. If the incom- 
mensupate term is absent [5 = 0), it is known from the exact Bethe ansatz solution and conformal field 
theorjO that the model in dl7j ) is gapless at zero magnetic field if — 1 < A < 1. If we now add a 
dimerization (q — and 5 is small), then the system becomes gapped and the gap scales as 

AE = 5 

where K = ^ AN . (18) 

2cos-!(-A) v ' 

[Actually, a gap opens up only if K < 2. If K > 2, the dimerization term is irrelevant in the sense of 
the renormalization group, and a gap is not generated]. It is useful to state the result in Eq. ( |l8| ) in 
the language of bosonization. For q — 0, the model in flyf ) is equivalent, at low energies, to a bosonic 
sine-Gordon theory described by the Lagrangian densityGT 



C = 

2ttK 



{dcj)/dt) 2 - (d(f>/dx) 2 1 - a S 2 / (2 - K ^ [ 1 - cos(20) ] , (19) 



5 



where a is a positive constant whose numerical value is not important here; the main point is to note 
the exponent of 5 in the coefficient o£-cos(2<jf>). (There are also factors of In 5 due to the presence of a 
marginal operator in the XX Z modelL3, but we will ignore such terms here). Thus, the incommensurate 
term which is linear in S in the original theory in ([l7]) becomes a cosine term with a different exponent for 
5 in the bosonic theory. This is because a nontriviaJ-cenormalization occurs in the process of deriving the 
low-energy bosonic theory from the fermionic oneo'til This renormalization occurs even if the fermions 
are noninteracting, i.e., for the XY spin-1/2 chain with A = in ( fl7| ) and K = 1; the reason for 
this is that the sine-Gordon theory is always strongly interacting. These strong interactions are also 
responsible for the large renormalizations of the correct quantum spectrum compared to the naive (i.e. 
classical) spectrum that one obtains from the sine-Gordon theory, namely, the classical soliton mass for the 
fermionic excitations and the quadratic fluctuations around <p — spectrum for the bosonic excitationso. 
The noninteracting fermionic model in ([!]) has no such renormalizations, which is why we did not use 
bosonization in the earlier part of this paper. 

Let us now change q slightly away from 0. Then bosonization again yields a theory of the sine-Gordon 
type, except that the coefficient of the cos(2</>) term in Eq. ( |l9| ) gets modified to produce 

- a \6sm{qx)\ 2n2 - K) [ 1 - cos{2(j>) ] , (20) 

where we have again shifted x to convert cos(<7x) to sin(qx). Unlike ([l9j), the theory in Eq. (^p| ) cannot be 
solved analytically, either in classical or in quantum mechanics. However we can make some qualitative 
statements about the low-energy spectrum if q/5 <C 1. When calculating the spectrum of small oscillations 
about a classical ground state <f> = 0, we find it convenient to scale x by a factor of 

a = (qd) 1 ^ 3 -^ . (21) 

This gives the eigenvalue equation 

_ g + AnKa \x\ 2 ^ 2 ~ K ^ = (^f c/> (22) 

for eigenstates lying close to the origin x = 0; we have approximated sin(qx) by qx and sin(20) by 2(f). 
Since ( p2| ) is the Schrodinger equation with a confining potential, we see that the energy can take several 
discrete values which are given by numerical factors multiplying a. [Note that our earlier results for 
noninteracting fermions agree with this scaling argument if we set K = 1] . These discrete values will 
then spread out into bands once we include tunneling between the different wells centered at the points 
x = nn. Similarly, we can find the energy of a classical soliton which goes from <f> = at x — > — oo to 
<fi = 7r at x — > oo; once again we can show by scaling that this will be given by a times some numerical 
factor. [Of course, all the numerical factors will get renormalized due to quantum corrections, but the 
power-law dependence on q5 is not expected to change] . We therefore see that all the low- lying excitations 
have energies of the order of (qS) 1 ^ 3 "^; this is much smaller than the gap which exists exactly at q = 
given by AE in Eq. (jl^), since we are assuming that q/S is small. We thus see that in the limit q — ► 0, 
there are many states which lie within the gap AE. 

We can use semiclassical arguments as in Eq. ([l3]) to estimate the density of states p(E) at energies 
much smaller than the q = gap. If E <C S 1 ^ 2 ^^, the wave function for that state lies in a region where 
we can approximate sin(qx) by qx. Ignoring various numerical factors, we then get 

P{E) = IE ^ V^ 2 - iW (2 ~ K) ©( E 2 - (Sqx) 2 ^ 2 -^ ) ~ ^ . (23) 

The specific heat at temperatures much lower than the q — gap therefore scales as T 3 ~ K /5, which is 
again much larger than the exponential dependence which occurs at q = 0. 



2ttK 



(d<j>/dt) 2 - {d<j>/dxf 
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Before ending this paper, we should remark that the peculiar difference between the limit q — > and 
q = only occurs near dimerization, i.e., near wave number tt. If we take the incommensurate term in 
Eq. (|l|) to be of the form S cos(Q + q)n, where Q/ir is equal to some simple rational number different 
from 1, then there is no such singularity at q — 0. Let us see this for the case of noninteracting fermions. 
Imagine filling up the Fermi sea up to a Fermi energy Ep = — cos kp where kp — Q/2; the Fermi velocity 
is then vp = sinfci?- We define the continuum Dirac field as 

tp(x) — i/jr{ x ) exp (ik,Fx) + iPl(x) exp (—ikpx) , (24) 
After dropping terms like exp(Hkpx) which vary rapidly, the Dirac equations take the form 
i ( d/dt + vpd/dx ) ip R + — exp (iqx — ikp) ipL = , 
i ( d/dt — vpd/dx ) ipi, + -z exp (—iqx + ikp) ip R ; = . (25) 



The solutions of this equation have the plane wave form 

(ip R (x,t)\ , p , / a R exp (ikx + iqx/2 ~ ikp/2) \ , . 

\i>L(x,t)J ~ exp 1 %m > \a L exp (ikx - iqx /2 + ik F /2) J ' [M) 



From this, we find that there is a single energy gap extending from E — (vpq — S)/2 to E = (vpq + 5)/2; 
the size of the gap is independent of q, assuming that q and S are both small. If S is held fixed and q is 
increased from zero, a state first appears at zero energy when q reaches the value 5/vp. It is instructive 
to express all this in the language of bosonization. Since K = 1, the sine-Gordon Lagrangian takes the 
form 

£ = — — (deb/dt) 2 - ^f(dcj)/dx) 2 - ^— [ 1 - cos(24> + qx) } , (27) 



where we have fixed thej-cp efficient of cos(2</>) in such a way that, for q = 0, the soliton mass, including 
the quantum corrections^, is exactly equal to the gap from zero energy, i.e., 5/2. It is convenient to scale 
x in the action (S — J dtdxC) by a factor of vp, and then shift (ft by qvpx/2. This gives 



C = 2, 



(deft/dt) 2 - (d(ft/dx - qvp/2) 2 



7rS 2 , , ,. , 

" ~w [ 1 ~ cos( 0) ] • 



(28) 



The Hamiltonian therefore contains a boundary term — (qvF /2ir)[(ft(oo) — cft(— oo)]; this is equal to —qvp/2 
in a one-soliton state. We thus see that if q is increased from zero, the energy of the one-soliton state 
becomes equal to that of the ground state (whose soliton number is zero) whenrflpj-gaches the value S/vp. 
[We note that similar results appear in the context of incommensurate crystals&Lnj] . 

Thus, if Q ^ 7r, there is a critical and nonzero value of q at which the gap closes at zero energy. This 
is very different from the earlier case with Q = n, where there is always a band of energies lying around 
zero energy, so that there is no gap at zero energy for any q ^ 0. (Our results disagree with Ref. 7 which 
argues that there is a nonzero critical value of q for Q = n). 

An interesting problem for future study may be to examine what happens if the incommensuration 
is strong, i.e., if 5 is comparable to or larger than 1. [Both the perturbative and the continuum Dirac 
theory (or bosonization) approaches would then fail]. It is possible that the pattern of bands and gaps 
would bepjam^much more complicated in that case, perhaps with a Devil's staircase or a point spectrum 
structureoau'cl. The nature of the wave functions, namely, whether they remain extended or become 
localized, would also be of interest. 
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Figure Captions 

1. The bands as a function of the energy E and the wave number q, both in units of <5, taking S = 0.05. 
The numbers 1, 3 and 5 labeling the three biggest gaps are explained in the text. 

2. A finer view of the bands as a function of E/S and q/S, for S = 0.05. A comparison with the 
second-order perturbation results is indicated by the six solid lines. 
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